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Abstract 

The calculation of the next-to-leading order SUSY-QCD correc- 
tions to the production of squarks, gluinos and gauginos at the Teva- 
tron is reviewed. The NL0 corrections stabilize the theoretical pre- 
dictions of the various production cross sections significantly and lead 
to sizeable enhancements of the most relevant cross sections for scales 
near the average mass of the produced massive particles. We dis- 
cuss the phenomenological consequences of the results on present and 
future experimental analyses. 
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1 Introduction 



The search for supersymmetric particles is among the most important en- 
deavors of present and future high energy physics. At the upgraded pp col- 
lider Tevatron, the searches for squarks and gluinos, as well as for the weakly 
interacting charginos and neutralinos, will cover a wide range of the MSSM 
parameter space 

The cross sections for the production of SUSY particles in hadron colli- 
sions have been calculated at the Born level already quite some time ago @. 
Only recently have the theoretical predictions been improved by calculations 
of the next-to-leading order SUSY-QCD corrections |3|-|5|]. The higher-order 
corrections in general increase the production cross section compared to the 
predictions at the Born level and thereby improve experimental mass bounds 
and exclusion limits. Moreover, by reducing the dependence of the cross sec- 
tion on spurious parameters, i. e. the renormalization and factorization scales, 
the cross sections in NLO are under much better theoretical control than the 
leading-order estimates. 

The paper is organized as follows. In Section 2 we shall review the calcu- 
lation of the next-to-leading order SUSY-QCD corrections [||-[5|], by using the 
case of qq production as an example. The NLO results for the production of 
squarks and gluinos are presented in Section 3. We first focus on the scalar 
partners of the five light quark flavors, which are assumed to be mass degen- 
erate. The discussion of final-state stop particles, with potentially large mass 
splitting and mixing effects, is presented in Section 4. In Section 5 we discuss 
the NLO cross sections for the production of charginos and neutralinos. We 
conclude the paper with a summary of the relevant MSSM particle produc- 
tion cross sections at the upgraded Tevatron, including next-to-leading order 
SUSY-QCD corrections.^] 

2 SUSY-QCD corrections 

The evaluation of the SUSY-QCD corrections consists of two pieces, the 
virtual corrections, generated by virtual particle exchanges, and the real cor- 
rections, which originate from real-gluon radiation as well as from processes 
with an additional massless (anti)quark in the final state. 



1 The MSSM Higgs sector will not be discussed here, see instead Ref.Q. 
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2.1 Virtual corrections 

The one-loop virtual corrections, i. e. the interference of the Born matrix ele- 
ment with the one-loop amplitudes, are built up by gluon, gluino, quark and 
squark exchange contributions (see Fig. [pa). We have adopted the fermion 
flow prescription for the calculation of matrix elements including Majo- 
rana particles. The evaluation of the one-loop contributions has been per- 
formed in dimensional regularization, leading to the extraction of ultraviolet, 
infrared and collinear singularities as poles in e = (4 — n)/2. For the chiral 
75 coupling we have used the naive scheme, which is well justified in the 
present analysis at the one-loop level.0 After summing all virtual correc- 
tions no quadratic divergences are left over, in accordance with the general 
property of supersymmetric theories. The renormalization of the ultravi- 
olet divergences has been performed by identifying the squark and gluino 
masses with their pole masses, and defining the strong coupling in the MS 
scheme including five light flavors in the corresponding (3 function. The mas- 
sive particles, i.e. squarks, gluinos and top quarks, have been decoupled by 
subtracting their contribution at vanishing momentum transfer ||. In di- 
mensional regularization, there is a mismatch between the gluonic degrees of 
freedom (d.o.f. = n — 2) and those of the gluino (d.o.f. = 2), so that SUSY 
is explicitly broken. In order to restore SUSY in the physical observables in 
the massless limit, an additional finite counter-term is required for the renor- 
malization of the novel qgq vertex. These counter-terms have been shown to 
render dimensional regularization consistent with supersymmetry ||. 

2.2 Real corrections 

The real corrections are generated by real-gluon radiation off all colored 
particles and by final states with an additional massless (anti) quark, obtained 
from interchanging the final state gluon with a light quark in the initial state 
(see Fig. |l|b). The phase-space integration of the final-state particles has been 
performed in n = 4 — 2e dimensions, leading to the extraction of infrared and 
collinear singularities as poles in e. After evaluating all angular integrals 
and adding the virtual and real corrections, the infrared singularities cancel. 
The left-over collinear singularities are universal and are absorbed in the 
renormalization of the parton densities at next-to-leading order. We have 
defined the parton densities in the conventional MS scheme including five 

2 We have explicitly checked that the results obtained with a consistent 75 scheme are 
identical to those obtained with the naive scheme. 
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light flavors, i.e. the squark, gluino and top quark contributions are not 
included in the mass factorization. We finally obtain an ultraviolet, infrared 
and collinear finite partonic cross section. 

There is, however, an additional class of physical singularities, which have 
to be regularized. In the second diagram of Fig. |l]b, an intermediate qg* state 
is produced, before the (off-shell) gluino splits into a qq pair. If the gluino 
mass is larger than the common squark mass, and the partonic cm. energy is 
larger than the sum of the squark and gluino masses, the intermediate gluino 
can be produced on its mass-shell. Thus the real corrections to qq production 
contain a contribution of qg production. The residue of this part corresponds 
to qg production with the subsequent gluino decay g — > qq, which is already 
contained in the leading order cross section of qg pair production, including 
all final-state cascade decays. This term has been subtracted in order to de- 
rive a well-defined production cross section. Analogous subtractions emerge 
in all reactions: if the gluino mass is larger than the squark mass, the con- 
tributions from g — > qq, qq have to be subtracted, and in the reverse case the 
contributions of squark decays into gluinos have to subtracted. 



3 Results 

In the following, we will present numerical results for SUSY particle pro- 
duction cross sections at the upgraded Tevatron (a/s = 2 TeV), including 
SUSY-QCD corrections. The hadronic cross sections are obtained from the 
partonic cross sections by convolution with the corresponding parton densi- 
ties. We have adopted the CTEQ4L/M parton densities |jlO| for the numerical 



results presented below. The uncertainty due to different parametrizations 
of the parton densities in NLO is less than ~ 15 %. The average final state 
particle mass is used as the central value of the renormalization and factor- 
ization scales and the top quark mass is set to m t = 175 GeV. The i^-factor 
is defined as K = unlo/o'lo, with all quantities parton densities, 

parton cross section) calculated consistently in lowest and in next-to-leading 
order. 



3.1 Production of Squarks and Gluinos 

Squarks and gluinos can be produced in different combinations via pp — > 
qq,qq,qg,gg- We first focus on the five light-flavored squarks, taken to be 
mass degenerate. At the central renormalization and factorization scale Q = 
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m, where m denotes the average mass of the final-state squarks/gluinos, 
the SUSY-QCD corrections are large and positive, increasing the total cross 
sections in general by 10-90% |J. This is shown in Fig. |2|, where the K- 
factors are presented as a function of the corresponding SUSY particle mass. 
The inclusion of SUSY-QCD corrections leads to an increase of the lower 
bounds on the squark and gluino masses by 10-30 GeV with respect to the 
leading-order analysis. 

The residual renormalization/factorization scale dependence in leading 
and next-to-leading order is presented in Fig. ^. The inclusion of the next- 
to-leading order corrections reduces the scale dependence by a factor 3-4 
relative to the lowest order and reaches a typical level of ~ 15%, when varying 
the scale from Q = 2m to Q = m/2. This may serve as an estimate of the 
remaining theoretical uncertainty due to uncalculated higher-order terms. 

Finally, we have evaluated the QCD-corrected single-particle exclusive 
transverse-momentum and rapidity distributions for all different processes. 
As can be inferred from Fig. [|, the modification of the normalized distribu- 
tions in next-to-leading compared to leading order is less than about 15% 
for the transverse-momentum distributions and even smaller for the rapid- 
ity distributions. It is thus a sufficient approximation to rescale the leading 
order distributions uniformly by the if-factors of the total cross sections. 

3.2 Stop Pair Production 

Stop production has to be considered separately since the strong Yukawa cou- 
pling between top/stop and Higgs fields gives rise to potentially large mixing 
effects and mass splitting. At leading-order in the strong coupling constant 
a s , only diagonal pairs of stop quarks can be produced in hadronic collisions, 
PP —> In contrast to the production of light-flavor squarks, the 

leading-order t-channel gluino exchange diagram is absent for stop production 
via qq initial states, since top quarks are not included in the parton densities. 
The leading-order stop cross section is thus in general significantly smaller 
than the leading order cross section for producing light-flavor squarks, where 
the threshold behavior is dominated by t-channel gluino exchange. Mixed 
tit 2 pair production can safely be neglected since it can proceed only via 
one-loop a s or tree level Gp amplitudes and is suppressed by several orders 
of magnitude The evaluation of the QCD corrections proceeds along 
the same lines as in the case of squarks and gluinos.0 The strong coupling 

3 The results obtained for the case of stop production can also be used to predict the 
sbottom pair cross section at NLO including mixing and mass splitting. 
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and the parton densities have been defined in the MS scheme with five light 
flavors contributing to their scale dependences, while the stop masses are 
renormalized on-shell. 

The magnitude of the SUSY-QCD corrections is illustrated by the K- 
factors at the central scale Q = mi in Fig. |5|. In the mass range relevant for 
the searches at the Tevatron, the SUSY-QCD corrections are positive and 
reach a level of 30 to 45% if the gg initial state dominates. If, in contrast, 
the qq initial state dominates, the corrections are small. The relatively large 
mass dependence of the i^-factor for stop production at the Tevatron can 
therefore be attributed to the fact that the gg initial state is important for 
small mi, whereas the qq initial state dominates for large m 4 ~. 

In complete analogy to the squark/gluino case, the scale dependence of 
the stop cross section is strongly reduced, to about 15% at next-to-leading 
order in the interval mi/ 2 < Q < 2m^. The virtual corrections at the NLO 
level depend on the stop mixing angle, the squark and gluino masses, and 
on the mass of the second stop particle. It turns out, however, that these 
dependences are very weak for canonical SUSY masses and can safely be 
neglected, as can be inferred from the light-stop production cross section 
in Fig. |5|. On the other hand, internal particles with masses smaller than 
the external particle mass, e.g. a light stop state propagating in the loops 
for heavy stop production, will contribute to the cross section. This feature 
explains the small but noticeable difference between the t\ and t<i i^-factors 
at mi = 300 GeV shown in Fig. |5|. 

The next-to-leading order transverse-momentum and rapidity distribu- 
tions are presented in Fig. || While the shape of the rapidity distribution is 
almost identical at leading and next-to-leading order, the transverse momen- 
tum carried away by hard gluon radiation in higher orders softens the NLO 
transverse momentum distribution considerably. 

3.3 Chargino and Neutralino Production 

At leading order, the production cross sections for chargino and neutralino 
final states depend on several MSSM parameters, i.e. Mi,M 2 ,n and tan/? 
[J2J. The cross sections are sizeable for chargino/neutralino masses below 
about 100 GeV at the upgraded Tevatron. Due to the strong sensitivity to 
the MSSM parameters, the extracted bounds on the chargino and neutralino 
masses depend on the specific region in the MSSM parameter space [|IJ. The 
outline of the determination of the QCD corrections is analogous to the 
previous cases of squarks, gluinos and stops. The resonance contributions 
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due to gq — > with g — > qxj have to be subtracted in order to avoid 
double counting of the associated production of electroweak gauginos and 
strongly interacting squarks. The parton densities have been defined with 
five light flavors contributing to their scale evolution in the MS scheme, while 
the t-channel squark masses have been renormalized on-shell.0 

At the average mass scale, the QCD corrections enhance the production 
cross sections of charginos and neutralinos typically by about 10-35% (see 
Fig. [7]), depending in detail on the final state and the choice of MSSM pa- 
rameters. The leading order scale dependence is reduced to about 10% at 
next-to- leading order (see Fig. ^), which implies a significant stabilization of 
the theoretical prediction for the production cross sections [|[. 

The individual leading order contributions of the s-channel gauge boson 
and the t, w-channel squark exchange are presented in Fig. For neutralino 
pair production the (t + -^-channel contributions are by far dominating, 
while the s-channel and interference terms are suppressed. Since the x? 2 
states are predominantly gaugino-like, this reflects the absence of a purely 
neutral trilinear gauge boson coupling in the Standard Model. Contrary to 
that the s-channel of xtx2 production is dominant and the (t + w)-channel 
term is suppressed for large squark masses. However, the interference turns 
out to be sizeable. 

4 Conclusions 

We have reviewed the status of SUSY particle production at the upgraded 
Tevatron at next-to-leading order in supersymmetric QCD. A collection of 
relevant sparticle production cross sections is shown in Fig. |9|. The higher- 
order corrections at the average mass scale of the massive final-state particles 
significantly increase the production cross section compared to the predic- 
tions at the Born level. Experimental mass bounds are therefore shifted 
upwards. Moreover, the theoretical uncertainties due to variation of renor- 
malization/factorization scales are strongly reduced to a level of typically 
~ 15 %, so that the cross sections in next-to-leading order SUSY-QCD are 
under much better theoretical control than the leading order estimates. The 
NLO results for total cross sections and differential distributions are available 

4 The next-to-leading order SUSY-QCD corrections to slepton pair production can 
be trivially obtained from the corresponding results for chargino/neutralino production. 
Numerically, the SUSY-QCD corrections for slepton production agree with the pure QCD 
corrections [ fj"T| , provided the squark and gluino mass are not chosen smaller than the final 
state slepton mass. 



7 



in the form of the computer code PROSPINO ||12|| . 
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Figure 1: Typical diagrams of the virtual (a) and real (b) corrections. 
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Figure 2: K-factors of the squark and gluino production cross sections at 
Q — rn. 
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Figure 3: Scale dependence of the total squark and gluino production cross 
sections for rriq = 250 GeV and m g = 300 GeV. 




Figure 4: Normalized transverse-momentum and rapidity distributions for 
squark production at Q = m. Mass parameters: rriq = 250 GeV and m g = 
300 GeV. 
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Figure 5: Left: K-factor of the stop production cross sections at Q = mi 
as a function of the stop masses [top/bottom scale]. Right: Production cross 
sections for the light stop state. The thickness of the NLO curve represents 
the dependence of the cross sections on the stop mixing angle and the gluino 
and squark masses. The shaded band indicates the theoretical uncertainty due 
to the scale dependence [m/2 < Q < 2m]. 
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Figure 6: Normalized transverse-momentum and rapidity distributions for 
stop production at Q = m t ~ = 200 GeV. 
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Figure 7: Left: K-factor of the X2X2 an d X1X2 production cross sections at 
the central scale Q = m. Right: Scale dependence of the X1X2 cross section. 
SUGRA parameters: m = 100 GeV, A = 300 GeV, tan/5 = 4, /z > 0. 
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Figure 8: The s-channel, (t + u)-channel and interference contributions to 
the leading order X2X2 an ^ X1X2 production cross section. SUGRA parame- 
ters as in Fig. [?| 
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Figure 9: The NLO production cross sections included in PROSPINO as a 
function of the final state particle mass; the arrows indicate the SUGRA 
inspired scenario: mi/ 2 = 150 GeV, m = 100 GeV, A = 300 GeV, tan/5 = 
4, /i > 0. All cross sections are given at the average mass scale of the massive 
final-state particles. 
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